The natural combination of full and image‐based waveform inversion          下载免费PDF全文

Tariq Alkhalifah  Zedong Wu 《Geophysical Prospecting》2016,64(1):19-30

Integrating migration velocity analysis and full waveform inversion can help reduce the high non‐linearity of the classic full waveform inversion objective function. The combination of inverting for the long and short wavelength components of the velocity model using a dual objective function that is sensitive to both components is still very expensive and have produced mixed results. We develop an approach that includes both components integrated to complement each other. We specifically utilize the image to generate reflections in our synthetic data only when the velocity model is not capable of producing such reflections. As a result, we get the migration velocity analysis working when we need it, and we mitigate its influence when the velocity model produces accurate reflections (possibly first for the low frequencies). This is achieved using a novel objective function that includes both objectives. Applications to a layered model and the Marmousi model demonstrate the main features of the approach.  相似文献   

Effective ellipsoidal models for wavefield extrapolation in tilted orthorhombic media     

Umair Bin Waheed  Tariq Alkhalifah 《Studia Geophysica et Geodaetica》2016,60(3):349-369

Wavefield computations using the ellipsoidally anisotropic extrapolation operator offer significant cost reduction compared to that for the orthorhombic case, especially when the symmetry planes are tilted and/or rotated. However, ellipsoidal anisotropy does not provide accurate wavefield representation or imaging for media of orthorhombic symmetry. Therefore, we propose the use of ‘effective ellipsoidally anisotropic’ models that correctly capture the kinematic behaviour of wavefields for tilted orthorhombic (TOR) media. We compute effective velocities for the ellipsoidally anisotropic medium using kinematic high-frequency representation of the TOR wavefield, obtained by solving the TOR eikonal equation. The effective model allows us to use the cheaper ellipsoidally anisotropic wave extrapolation operators. Although the effective models are obtained by kinematic matching using high-frequency asymptotic, the resulting wavefield contains most of the critical wavefield components, including frequency dependency and caustics, if present, with reasonable accuracy. The proposed methodology offers a much better cost versus accuracy trade-off for wavefield computations in TOR media, particularly for media of low to moderate anisotropic strength. Furthermore, the computed wavefield solution is free from shear-wave artefacts as opposed to the conventional finite-difference based TOR wave extrapolation scheme. We demonstrate applicability and usefulness of our formulation through numerical tests on synthetic TOR models.  相似文献   

Spatial Periodic Boundary Condition for MODFLOW     

Tariq Laattoe  Vincent E.A. Post  Adrian D. Werner 《Ground water》2014,52(4):606-612

Small‐scale hyporheic zone (HZ) models often use a spatial periodic boundary (SPB) pair to simulate an infinite repetition of bedforms. SPB's are common features of commercially available multiphysics modeling packages. MODFLOW's lack of this boundary type has precluded it from being effectively utilized in this area of HZ research. We present a method to implement the SPB in MODFLOW by development of the appropriate block‐centered finite‐difference expressions. The implementation is analogous to MODFLOW's general head boundary package. The difference is that the terms on the right hand side of the solution equations must be updated with each iteration. Consequently, models that implement the SPB converge best with solvers that perform both inner and outer iterations. The correct functioning of the SPB condition in MODFLOW is verified by two examples. This boundary condition allows users to build HZ‐bedform models in MODFLOW, facilitating further research using related codes such as MT3DMS and PHT3D.  相似文献   

Uncertainty in recurrence rates of large magnitude events due to short historic catalogs   总被引：2，自引：0，他引：2  

Tariq Mohammed  Gail M. Atkinson  Karen Assatourians 《Journal of Seismology》2014,18(3):565-573

Seismic hazard analysis requires knowledge of the recurrence rates of large magnitude earthquakes that drive the hazard at low probabilities of interest for seismic design. Earthquake recurrence is usually determined through studies of the historic earthquake catalogue for a given region. Reliable historic catalogues generally span time periods of 100–200 years in North America, while large magnitude events (M?≥?7) have recurrence rates on the order of hundreds or thousands of years in many areas, resulting in large uncertainty in recurrence rates for large events. Using Monte Carlo techniques and assuming typical recurrence parameters, we simulate earthquake catalogues that span long periods of time. We then split these catalogues into smaller catalogues spanning 100–200 years that mimic the length of historic catalogues. For each of these simulated “historic” catalogues, a recurrence rate for large magnitude events is determined. By comparing recurrence rates from one historic-length catalogue to another, we quantify the uncertainty associated with determining recurrence rates from short historic catalogues. The use of simulations to explore the uncertainty (rather than analytical solutions) allows us flexibility to consider issues such as the relative contributions of aleatory versus epistemic uncertainty, and the influence of fitting method, as well as lending insight into extreme-event statistics. The uncertainty in recurrence rates of large (M?>?7) events is about a factor of two in regions of high seismicity, due to the shortness of historic catalogues. This uncertainty increases greatly with decreasing seismic activity. Uncertainty is dependent on the length of the catalogue as well as the fitting method used (least squares vs. maximum likelihood). Examination of 90th percentile recurrence rates reveals that epistemic uncertainty in the true parameters may cause recurrence rates determined from historic catalogues to be uncertain by a factor greater than 50.  相似文献   

Sea Level Variations and Geomorphological Changes in the Coastal Belt of Pakistan     

Tariq Masood Ali Khan  D. A. Razzaq  Qamar-Uz-Zaman Chaudhry  Dewan Abdul Quadir  Anwarul Kabir  Majajul Alam Sarker 《Marine Geodesy》2013,36(1-2):159-174

The UNEP in its regional seas program in 1989 has included Pakistan in a group of countries which are vulnerable to the impact of rising sea level. If the present trend of sea level rise (SLR) at Karachi continues, in the next 50 years the sea level rise along the Pakistan Coast will be 50 mm (5 cm). Since the rising rates of sea level at Karachi are within the global range of 1-2 mm/year, the trends may be treated as eustatic SLR. Historical air temperature and sea surface temperature (SST) data of Karachi also show an increasing pattern and an increasing trend of about 0.67°C has been registered in the air temperature over the last 35 years, whereas the mean SST in the coastal waters of Karachi has also registered an increasing trend of about 0.3°C in a decade. Sindh coastal zone is more vulnerable to sea level rise than Baluchistan coast, as uplifting of the coast by about 1-2 mm/year due to subduction of Indian Ocean plate is a characteristic of Baluchistan coast. Within the Indus deltaic creek system, the area nearby Karachi is more vulnerable to coastal erosion and accretion than the other deltaic region, mainly due to human activities together with natural phenomena such as wave action, strong tidal currents, and rise in sea level. Therefore, The present article deals mainly with the study of dynamical processes such as erosion and accretion associated with sea level variations along the Karachi coast and surrounding Indus deltaic coastline. The probable beach erosion in a decade along the sandy beaches of Karachi has been estimated. The estimates show that 1.1 mm/year rise in sea level causes a horizontal beach loss of 110 mm per year. Therefore, coast eroded with rise in sea level at Karachi and surrounding sandy beaches would be 1.1 m during a period of next 10 years. The northwestern part of Indus delta, especially the Gizri and Phitti creeks and surrounding islands, are most unstable. Historical satellite images are used to analyze the complex pattern of sediment movements, the change in shape of coastline, and associated erosion and accretion patterns in Bundal and Buddo Islands. The significant changes in land erosion and accretion areas at Bundal and Buddo Islands are evident and appear prominently in the images. A very high rate of accretion of sediments in the northwestern part of Buddo Island has been noticed. In the southwest monsoon season the wave breaking direction in both these islands is such that the movement of littoral drift is towards west. Erosion is also taking place in the northeastern and southern part of Bundal Island. The erosion in the south is probably due to strong wave activities and in the northeast is due to strong tidal currents and seawater intrusion. Accretion takes place at the northwest and western parts of Bundal Island. By using the slope of Indus delta, sea encroachment and the land area inundation with rising sea level of 1 m and 2 m have also been estimated.  相似文献   

Investigation of the groundwater modelling component of the Integrated Water Flow Model (IWFM)     

Ali Ercan  Emin C. Dogrul  Tariq N. Kadir 《水文科学杂志》2013,58(16):2834-2848

ABSTRACT

The Integrated Water Flow Model (IWFM), developed by the California Department of Water Resources, is an integrated hydrological model that simulates key flow processes including groundwater flows, streamflow, stream–aquifer interactions, rainfall–runoff and infiltration. It also simulates the agricultural water demand as a function of soil, crop and climatic characteristics, as well as irrigation practices, and allows the user to meet these demands through pumping and stream diversions. This study investigates the modelling performance of the groundwater module of IWFM using several hypothetical test problems that cover a wide range of settings and boundary conditions, by comparing the simulation results with analytical solutions, field and laboratory observations, or with results from MODFLOW outputs. The comparisons demonstrate that IWFM is capable of simulating various hydrological processes reliably.